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Abstract 



This paper establishes a blowup criterion for the three-dimensional viscous, 
compressible, and heat conducting magnetohydrodynamic (MHD) flows. It is es- 
sentially shown that for the Cauchy problem and the initial-boundary-value one 
of the three-dimensional compressible MHD flows with initial density allowed to 
vanish, the strong or smooth solution exists globally if the density is bounded from 
above and the velocity satisfies the Serrin's condition. Therefore, if the Serrin norm 
of the velocity remains bounded, it is not possible for other kinds of singularities 
(such as vacuum states vanish or vacuum appears in the non-vacuum region or 
even milder singularities) to form before the density becomes unbounded. This 
criterion is analogous to the well-known Serrin's blowup criterion for the three- 
dimensional incompressible Navier-Stokes equations, in particular, it is indepen- 
dent of the temperature and magnetic field and is just the same as that of the 
barotropic compressible Navier-Stokes equations. As a direct application, it is 
shown that the same result also holds for the strong or smooth solutions to the 
three-dimensional full compressible Navier-Stokes system describing the motion of 
a viscous, compressible, and heat conducting fluid. 

Keywords: compressible magnetohydrodynamic system, full compressible Navier- 
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1 Introduction 

In this paper, we consider the system of partial differential equations for the three- 
dimensional viscous, compressible, and heat conducting magnetohydrodynamic (MHD) 
flows in the Eulerian coordinates |21j 

/■ 

pt + dw{pu) = 0, 
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{pu)t + div(/ju (g) u) 
c^ipO) t + div {puO)] 



pAu -{p + A)Vdivu + VP = (curl H) x H, 
kAO + Pdwu = 2^|D('u)p + A(div?i)2 + u\cutI H 
'AH, divH = 0, 
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where t > is time, a; G C is the spatial coordinate, and p,u = {ui,U2,u^)^^ , 
0, P = Rp9 {R > 0), and H = (Hi, H2, H3) , represent respectively the fluid density, 
velocity, absolute temperature, pressure, and magnetic field; is the deformation 

tensor given by 

D(n) = l(Vn + (V^xn. 

The constant viscosity coefficients fi and A satisfy the physical restrictions 

/i>0, 2/i + 3A>0. (1.2) 

Positive constants c^,, k, and ly are respectively the heat capacity, the ratio of the heat 
conductivity coefficient over the heat capacity, and the magnetic diffusivity acting as a 
magnetic diffusion coefficient of the magnetic field. 

The equations (jl.ip will be studied with initial condition: 

{p,u,0,H){x,O) = {po,uo,0Q,Ho)ix), X G 0, (1.3) 

and one of the following boundary conditions: 

1) If $7 = R^, for constant p > 0, {p,u,6, H) satisfies the far field condition: 

{p,u, H,6){x,t) ^ {p, 0,0,0) as |x| — )- 00; (1.4) 

2) If is a bounded smooth domain in M'^, (n, 6, H) satisfies 

u = 0, T- = 0, H = ondn, (1.5) 
on 

where n = (ni,n2,?T.3) is the unit outward normal to 90. 

The compressible MHD system (jl.ip is a combination of the compressible Navier- 
Stokes equations of fluid dynamics and Maxwells equations of electromagnetism. In- 
deed, the equations (jl.ip i . (|l.ip 9. and (jl.iP ':; describe, respectively, the conservation of 
mass, momentum, and energy. In addition, it is well-known that the electromagnetic 
fields are governed by Maxwells equations. In magnetohydrodynamics, the displace- 
ment current can be neglected ( [2T]). As a consequence, the equation (jl.ip zi is called 
the induction equation, and the electric field can be written in terms of the magnetic 
field H and the velocity u, 

E = iyVxH-uxH. 

Although the electric field E does not appear in the compressible MHD system (|l.ip , it 
is indeed induced according to the above relation by the moving conductive flow in the 
magnetic fleld. In particular, when there is no electro-magnetic effect, that is, H = 0, 
the compressible MHD system (jl.ip reduces to the following full compressible Navier- 
Stokes system describing the motion of a viscous, compressible, and heat conducting 
fluid: 

pt + div(p'u) = 0, 

< {pu)t + div(pM (g) u) - pAu -{p + A)Vdivu + VP = 0, (1.6) 
^c,a[{p9)t + divipuO)] - kA6 + Pdivu = 2/i|D(M)|2 + A(divu)2. 

There is a considerable body of literature on the multi-dimensional full compressible 
Navier-Stokes system ()1.6p and compressible MHD one (jl.ip by physicists and mathe- 
maticians because of their physical importance, complexity, rich phenomena, and math- 
ematical challenges; see (ilSElElIinillllliainiEaiHEaETlIMllMlET^ and the references 
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cited therein. However, many physically important and mathematically fundamental 
problems are still open due to the lack of smoothing mechanism and the strong non- 
linearity. For example, although the local strong solutions to the compressible MHD 
system (jl.ip with large initial data were respectively obtained by [36J and |7J in the 
cases that the initial density is strictly positive and that the density is allowed to vanish 
initially, whether the unique local strong solution can exist globally is an outstanding 
challenging open problem. 

Therefore, it is important to study the mechanism of blowup and structure of possible 
singularities of strong (or smooth) solutions to the compressible MHD system (jl.ip and 
to the full compressible Navier-Stokes one (jl.6p . The pioneering work can be traced 
to Serrin's criterion [29j on the Leray-Hopf weak solutions to the three-dimensional 
incompressible Navier-Stokes equations, which can be stated that if a weak solution u 
satisfies 

u G L'(0, T;U), - + -<!, 3 < r < oo, (1.7) 
s r 

then it is regular. Later, He-Xin [9] showed that the Serrin's criterion ()1.7p still holds 
even for the strong solution to the incompressible MHD equations. 

Recently, Huang-Li-Xin [T7j extended the Serrin's criterion (|1.7p to the barotropic 
compressible Navier-Stokes equations and showed that if T* < oo is the maximal time 
of existence of a strong (or classical) solution {p,u), then 

^lim^ (||divn||ii(o,T;L-) + ll^^llL=(o,r;L'-)) = oo, (1.8) 

and 

^1™. (IIpIIl«=(0,T;L°°) + \\u\\l''{0,T;L-)) = OO, (1.9) 

with r and s as in p.7p . For more information on the blowup criteria of barotropic 
compressible flow, we refer to [12l[T3l[I71[l8l[20l[3l] and the references therein. Later 
Xu-Zhang [38] extended the results of [l7j to the isentropic compressible MHD system 
and obtained that the same blow-up criterion (ll.9p holds. 

When it comes to the full compressible Navier-Stokes system (II. 6p . the problem is 
much more complicated. Let T* < oo be the maximal time of existence of a strong (or 
classical) solution {p,u,9) to the system (II. 6p . Besides (jl.2p . under the condition that 

7p > A, (1.10) 

Fan-Jiang-Ou [6] obtained that 

ylim,(l|6'||L-=(0,T;L-=) + I|Vu||l1(o,T;L°°)) = OO. 

Recently, under just the physical restrictions (II. 2p . Huang-Li [14] and Huang-Li-Xin [18] 
established the following blowup criterion: 

^hm^ (||6'||l2{o,T;L-) + I|2)('")IIli(o,T;L-)) = oo, 

where D{u) is the deformation tensor. Later, in the absence of vacuum. Sun- Wang- 
Zhang [35] showed that 

^lim. (||^IIl-(0,T;L-) + II {P,P'^) IL-(0,T;L-)) = 
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provided that (ll.2p and (jl.lOp both hold. Very recently, under just the physical restric- 
tions (jl.2p and allowing the initial density to vanish, Huang-Li- Wang |16j improved all 
the previous results [6|[n i [T8 | [35] by obtaining that (II. 8p still holds. It should be noted 
here that (jl.9p is much stronger than (jl.Sp and that whether (jl.Op holds or not remains 
open. 

For the compressible MHD system p.ip . let T* < oo be the maximal time of existence 
of a strong (or classical) solution {p,u,9, H). Lu et al [23j obtained that 



and 



^1™:,JIIpIIl°°(0,T;L°°) + II^IIl°°{0,T;L°°) + II Vn||i4(o,T;L2)) = OO, 



^hmJ||divu||ioc(o,T;L-=) + I|6'||l°°(o,T;L->) + l|Vu||i4(o,r;L2)) = oo. 



while Chen-Liu [S] showed that 

^lim J||Vu||/,i(o,T;L°°) + II^IIl°°(0,T;L°°)) = C»- 

The aim of this paper is to improve all the previous blowup criterion results on 
both the compressible MHD system (II. ip and the full compressible Navier-Stokes one 
(|1.6p by allowing initial vacuum states, and by describing the blowup mechanism just 
in terms of the Serrin-type criterion, (II. 9p . Before stating our main result, we first 
explain the notations and conventions used throughout this paper. We denote 

fdx = / fdx. 
Jn 

For 1 < p < oo and integer A; > 0, the standard homogeneous and inhomogeneous 
Sobolev spaces are denoted by: 

LP = LP{n), W'^^P = W'''P{n), D^^P = {-u G L\^^{n) I V^u e LP} , 
Dl= {u£ L^\Vu€ L^,u = on dn} , = L^ n D^, H'' = W^''^, 

^2,2 _ / e ^^'^ n 1)2,2 I V6I • n = on dO.) , for bounded f], 
\z)ini?2.2, forJ7 = IR3. 

Then, the strong solutions to the initial-boundary-value problem (|l.ip - (|1.3p together 
with ()1.4p or ()1.5p are defined as follows. 

Definition 1.1 (Strong Solutions). For p>0 and 9 = 0, {p,u,9, H) is called a strong 
solution to (jl.ip in Q X (0, T), if for some qq > 3, 

P>o, p-pe C([o,r];/7inw^i''?n), pt & C{[o,t];L^ n l"'^), 

{H, u) E ^([0, T];Dl n D^^^) n L'^{0, T; D^^'^o), H G C([0, T];H^) 
9>0, 9e C([0,r];Z)2;2) nL2(0,r;Z)2.9o), 

{Ht,ut,9t) G L2(0,r;Z)i'2), {Ht, ^m, ^9t) G L^{0,T; L^), 

and {p, u, 9, H) satisfies both (jl.ip almost everywhere in Q x (0, T) and ()1.3p almost 
everywhere in Q. 



Our main result can be stated as follows: 
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Theorem 1.1. For q G (3,6], assume that the initial data (po ^ O,tto,^o ^ 0, i^o) 
satisfies 



Po\uo\^ + PoO^o L\ HoGH^nH^ divi^o = 0, 
and the compatibility conditions 



(1.11) 



- pAuo -{p + A)Vdiv^io + RV{po9o) - (cmlHo) x Ho = Vpo5i, (1-12) 
kA9o + ^|Vno + (Vno)*n^ + A(divuo)2 + z^|curlFop = VA)52, (1-13) 

mf/i G L^- Let {p,u,9, H) be the strong solution to the initial boundary value 

problem (jl.ip - (jl.3p together with (jl.4p or (jl.Sp . If T* < oo is the maximal time of 
existence, then for r and s as in ()1.7|) , 

ylim,(llp||L-'(o,r;L°°) + II'"IIl''(o,T;L'-)) = oo- (1-14) 

If H = Hq = 0, Theorem 11.11 directly yields the following Serrin-type blowup criterion 
for the three-dimensional full compressible Navier-Stokes system (jl.6p . 



Theorem 1.2. For constants q G (3,6] and p > 0, assume that {po > 0, no,0o ^ 0) 
satisfies 

Po-p£H^nW^'^, uo^DlnD^'^, ^oG^o,n, po\uo\'' + poO^ e L\ 
and the compatibility conditions 

-pAuo -{p + A)Vdivno + RV{pq9q) = ^/p^gi, 

kA9o + ||Vno + (Vuo)*n^ + A(divuo)^ = Va)52, 

with gi,g2 G L^- Let {p,u,9) be the strong solution to the full compressible Navier- 
Stokes system p.6p together with 

{p, u, 9){x, 0) = {po, uo, 9o), X eQ, (1.15) 

and either for = M^, 

{p, u, 9) — )• (p, 0, 0) as \x\ — )• cxD, (1-16) 
or for a bounded smooth domain C M^, 

M = 0, — = on 50. (1.17) 

on 

IfT* < oo is the maximal time of existence, then 

^1™.(IIpIIl°°{0,T;L°°) + II^^IIl-(0,T;L'-)) = OO, (1-18) 

with r and s as in (jl.Zp . 
A few remarks are in order: 

Remark 1.1. The conclusion in Theorem \l.l\ is somewhat surprising since the criterion 
()1.14p is independent of the temperature and magnetic fields and just the same as those 
of barotropic compressible Navier-Stokes equations ( Jj7[ /J. 
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Remark 1.2. In 116, Theorem 1], we obtained that (II. Sh holds for the Cauchy problem 
of the full compressible Navier-Stokes system (jl.6p . Thus, 

lim ||divM||ii(o,T;L-') = oo, (1.19) 

provided that 

sup \\u\\LrmT;L=) < OO, (1.20) 

o<r<r* 

for r, s as in ()1.7p . It follows from the continuity equation ()1.6P i that for t G [0,T*), 

p{x,t) = poiy{0;x,t))exp i^- j divu{y{s;x,t),s)ds^ , (1.21) 
where y{s;x,t) is the characteristic curve defined by 

-j-y = u{y^s), y{t;x,t)=x. 
as 

The combination of p.l9p with p.2ip implies that there may hold for the density: 

1) The density remains bounded, that is, 

^1™. II^IIl°°(o,T;L-) < oo; (1-22) 

2) The density may concentrate, that is, 

^lim^ \\p\\l^{q,T;L^) = oo; (1.23) 

3) Vacuum states may vanish: There exists some xi £ 0, and xi{t) satisfying 
Poixi) = and y(0; xi{t),t) = x\ such that 

^Hm p(xi(t),t) > Co > 0; (1.24) 

J^) Vacuum states may appear in the non-vacuum region: There exists some X2 G ^ 
and X2{t) satisfying po{x2) > and y{0; X2{t),t) = X2 such that 

^lim p(x2(t),t) = 0. (1.25) 



Then one may ask: Which one or some of (jl.22p -~ (jl.25p will happen ? Theorem [77 
gives an answer to this question by obtaining that the density will concentrate provided 
that (|1.2Up holds. In other words, if the Serrin norm of the velocity remains bounded, 
it is not possible for other kinds of singularities (such as vacuum states vanish or 
vacuum appears in the non-vacuum region or even milder singularities) to form before 
the density becomes unbounded. Moreover, (jl.lSp still holds for the initial-boundary- 
value problem (jl.6p (jl.l5p (jl.l7p . Thus, Theorem greatly improves all the previous 
blowup criterion for the full compressible Navier-Stokes system ()1.6p 14 \1(}1 \35 ^ . 

Remark 1.3. If Q is a bounded smooth domain o/M^, Theorems \1.1\ and \1.2\ still hold 
when the boundary condition VO ■ n\gQ = is replaced by 6\qq = 0. 

Remark 1.4. Theorems \l.l\ and \1.2\ also hold respectively for classical solutions to the 
three-dimensional compressible MHD system (jl.ip and to the full compressible Navier- 
Stokes one (11.61). 
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We now comment on the analysis of this paper. 

Let {p,u,9,H) be a strong solution described in Theorem ll.il Suppose that ()1.14p 
were false, that is, 

^lim, {\\p\\l^(o,T;L°°) + \\u\\L''io,T;Ln) < M) < +00. (1.26) 

We want to show that 

sup (11/9 - /5||//invi/i>9 + ll^^lli/i + ll^^ll-ffi + ll-f^ll/fO <C< +00. 

0<t<T* 

Since the methods in all previous works [3l[6l[T6 t [23 t [35 t [55] depend crucially on either 
the L^L!^-norm of the temperature 9 or the LjL^-norm of the divergence of the 
velocity divn, some new ideas are needed to recover all the a priori estimates just under 
the assumption ()1.26p without any a priori bounds on the temperature, the magnetic 
field, and the divergence of the velocity. In fact, we prove (see Lemma [3. 3 p that a control 
of the Serrin norm of the velocity and L^L^-norm of the density implies a control on 
the L'^L^ norm of Vu. In order to obtain this control, the key observation is that, 
instead of the temperature 6, we treat the total energy E = c^O + which in turn 

greatly reduces the difficulties arising from the high nonlinearities of the temperature 
equation, (jl.ip q. Indeed, multiplying the equation of the conservative form of the total 
energy E (see (13. Sp ) by E yields that to bound the L^L^-norm of VE (see (13. 4p ). it is 
enough to control that of |?i||Vu|, which in fact can be reduced to the estimate of the 
L^L^-norm of Vu (see (|3.30p ). Then, to overcome the difficulty caused by the boundary 
when n is bounded, motivated by [ISIM], we decompose the velocity into two parts 
(see (j3.14p and (j3.18p ) which together with the L^'-estimate for the Lame system yield 
the desired bound on the L^L^-norm of Vn (see ()3.3ip ). Finally, the a priori estimates 
on both the L^L^-norm of the density gradient and the L\L'^-noim of the velocity 
gradient can be obtained simultaneously by solving a logarithm Gronwall inequality 
based on a logarithm estimate for the Lame system (see Lemma l2.3p and the a priori 
estimates we have just derived. 

The rest of the paper is organized as follows: In the next section, we collect some 
elementary facts and inequalities that will be needed later. The main result. Theorem 
ll.lt is proved in Section [3l 

2 Preliminaries 

In this section, we recall some known facts and elementary inequalities that will be 
used later. 

First, the following existence and uniqueness of local strong solutions when the initial 
density may not be positive and may vanish in an open set can be proved in a similar 
way as in [4] (cf. [7J). 

Lemma 2.1. Assume that the initial data {pQ > 0,uq,6q > 0,Hq) satisfy (jl.lip - (|1.13p . 
Then there exists a positive time Ti G (0, oo) and a unique strong solution (p, u, 9, H) to 
the initial-boundary-value problem p.ip - ()1.3p together with ()1.4p or ()1.5p on x (0, Ti]. 

Next, the following well-known Sobolev inequality will be used later frequently (see 
[26]). 
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Lemma 2.2. For p £ (l,oo) and q G (3,cxd), there exists a generic constant C > 0, 
which depends only on p, q such that for f G Dq and g £ Ci D^''^, we have 



1^6 <C||V/||i2, Ml^ <C\\g\\LP + C\\Vg\\L.. (2.1) 

Finally, we consider the following Lame system 

- nAv{x) - + X)VdiYv{x) = f{x), x e n, (2.2) 

where v = {vi,V2,V3), f = (/i,/2,/3), and fj,,X satisfy ()1.2p . The system ()2.2p is 
imposed on one of the following boundary conditions: 

1) Cauchy problem: $7 = M^, and 

v{x) — )• 0, as |x| — )• oo; (2.3) 

2) Dirichlet problem: O is a bounded smooth domain in M^, and 

V = on dn. (2.4) 

The following logarithm estimate for the Lame system (12.20 will be used to estimate 
||Vu||loo and WVpHi^nLi- 

Lemma 2.3. Let fi,X satisfy (II. 2p . Assume that f = divg where g = {gkj)3x3 with 
gkj G LP' r\ U r\ D^''^ for k,j = !,••• ,3, r G (l,oo), and q G (3, oo). Then the Lame 
system (12.20 together with (12.30 or (12.40 has a unique solution v G Dq n D^'^ n D^''^, 
and there exists a generic positive constant C depending only on fj., A, q, and r (besides 
Vt when Vt is hounded) such that 

\\yv\\Lr <C\\g\\Lr, (2.5) 

and 

\Wv\\loo <C{1 + ln(e + IIV^IIlOIIsIIl- + II^IIlO • (2-6) 
Proof. First, if = M^, direct calculations show that v = {vi,V2,V'^) with 

is the unique solution to the Cauchy problem (|2.2p (|2.3p and satisfies (|2.5p . 

Then, if is a bounded smooth domain of M^, it follows from |33j that the Dirichlet 
problem ()2.2p ()2.4p is of Petrovsky type. In Petrovsky's systems, roughly speaking, 
different equations and unknowns have the same "differentiability order", see [321 p. 126]. 
We also recall that Petrovsky's systems are an important subclass of Agmon-Douglis- 
Nirenberg (ADN) elliptic systems( |l]), having the same good properties of self-adjoint 
ADN systems. It follows from Solonnikov [321 Theorem 1.1] and |331 Theorem 5.1] that 
the solution v to the system ()2.2p together with (12. 4p can be represented as 

Vi{x) = j Gij{x,y)fj{y)dy, for ah x G fi, (2.7) 

by means of the Green function Gij = Gij{x,y) G C°°(Q x Q, \ D) with D = {{x,y) G 
X n|x = y} which satisfies that for every multi-indexes a = (01,02,03) and /3 = 
(/5i, /32) /^s) there is a constant Ca,p such that for all {x, y) G QxQ,\D, and i,j = l,--- ,3, 

\d^d^Gij{x,y)\<C^,p\x-y\-'-^"^-^('\ (2.8) 
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where |q| = 01+02 + 03 and /3 = /3i + /32 + /33. Moreover, the estimate (12. 5p is standard. 

Finally, it remains to prove (12. 6p . We will only deal with the Dirichlet problem (I2.2p 
()2.4p . since the same procedure holds for the Cauchy problem (j2.2p ()2.3p . Motivated 
by Beale-Kato-Majda [2], we introduce a small parameter 6 G (0, 1] which depends on 
V, and which will be fixed later. Using 6, we define a cut-off function %(s) satisfying 
7/5(5) = 1 for < s < (5, 775(5) = for s > 25, and \'ns'\s)\ < C5^^. It thus follows from 
(1231) that 



i'i(x) = / {r]5{\x-y\) + {l-r}s{\x-y\)))Gij{x,y)dkgkj{y)dy 

m{\x-y\)Gij{x,y)dkgkj{y)dy + j dy^r]s{\x - y\)Gij{x,y)gkjiy)dy 
(1 - r]s{\x - y\))dy^Gij{x,y)gkjiy)dy, 



where in the second equality we have used integrations by parts due to the fact that 
Gij{x,y)\QQ = for each x £ Hence, we have 

\Vv{x)\ <C I {Ws\\G,,\+r]s\V^Gij\) \Vg\dy 

+ G j {\4\\G^J\ + |??^||V,.G,,| + Ws\\VyG,j\) \g\dy (2.9) 

+ G j {l-m)\V^VyG^j\\g\dV- 
Each term on the right-hand side of (j2.9p can be estimated by (I2.8P as follows: 

W5\\Gij\+m\^xGij\) \Vg\dy 



and 



<C f (rV - 7/|-i + |x - y|-2) \Vg\dy 

J^r\{y\\x-y\<28} (2.10) 
/ 1-25 1-25 \ C?-!)/? 

< C / s-«/(«-l)s2rf5+ / IIV^IIl. 

(k'HG.,! + |r?^||V,G,,| + |7?^||V,G,,|) \g\dy 

j-25 

<GJ {6~^s~^ + 6-h''^)s'^ds\\g\\Loo 
< G\\g\\L^, 

(1 - V5)\'^x'^yGij\\g\dy 



(2.11) 



<C / + / \x - y\-^\giy)\dy 

\Jnn{y\s<\x-y\<i} Jnn{y\\x-y\>i}J (2.12) 

<C/ fi-^^s^dsll^llioo +C7 ( / s-^'/^'-^h^ds) \\g\\Lr 
< -C\\g\\L^ ln6 + G\\g\\Lr. 
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It follows from (|2:9ll- (|2l^ that 

\\Vv\\l^ < C (s^^'^^/^VgU. + (1 - lnS)\\g\\L^ + MlA . (2.13) 



Set 6 = mill 1 1, llVgW^f^'^ • Then (I2J3I1 becomes <^M)- We finish the proof of 
Lemma 12.31 

3 Proof of Theorem 11.11 

Before proving Theorem 11.11 we state some a priori estimates under the condition 
(ll.26p . First, we have 

Lemma 3.1. Under the condition (ll.26p . it holds that for q G [2, 12] and < T < T* , 

\\H\\lo^(o,T;Li) + £ j iHl^-^lVHl^dx < C, (3.1) 

where (and in what follows) C and Ci{i = 1, ■ ■ ■ ,6) denote generic constants depending 
only on Mq, fi, A, R,k, c^, T*,v, and the initial data (besides for hounded Q). 

Proof. Similar to [38j, multiplying (1.1)4 by q\H\'^^^H and integrating the resulting 
equation over yield that 

-| J \H\'idx + i^ J {q\H\'i-^\VH\'^ + q{q - 2)\H\'i-^\V\H\\'^) dx 

= - j q\H\'i~^ (^H-VH-u- • V|/7p^ dx 

_q{q_2)_ I \H\'i~'i^H-V\H\'){u-H)da 

< q\H\'^-^\VHfdx + Cq^ J \u\^\H\ 



(3.2) 



I'^dx 



< q\Hr^\VH\^dx + C\\u\\lr\\\H\'^/^\f,^r'^/^ 

q\H\''~'^\VHfdx + C6\\V\H\'^/'^\\l2+C{6){l + \\u\\ir)\\H\ 



Li- 



Choosing S suitably small in (13. 2p . we obtain (13. Ih directly after using Gronwall's 
inequality and (jl.26p . We thus finish the proof of Lemma 13.11 

Then, we derive the following key estimate on the specific energy E defined by 

E^cJ + ^-. (3.3) 
Lemma 3.2. Under the condition (|1.26p . it holds that 

j pE'^dx + K\\VE\\l2 <C j \u\^ [pE"^ + \Vu\^)dx + Ci\\VH\\]jr 



+ C\\Vu\\l2+C / pE'^dx. 



(3.4) 
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pE'^dx + K / \VE\'^dx 



Proof. First, it follows from (II. ip that E satisfies 

{pE)t + div{pEu) - —AE = divF - WWdiU^ + -ji^l^divu + u\cm\H\^, (3.5) 
Cy 2 

with 

F ^ |^V(|np) + pu-Vu + Audivu -Pu + {u- H)H - ]^\H\'^u. 

Next, applying standard maximum principle to (jl.iP o together with ^ (c.f. [6l[8]) 
shows 

inf (9(x,t)>0. 

K3x[0,T] 

Multiplying (j3.5p by c^i? and integrating the resulting equality over fi, we obtain after 
integration by parts and using (jl.ip -^ that 

Cy d 

<C j {\u\\Vu\+ pe\u\)\VE\dx + C j {\u\\H\^\VE\ + \Vu\\H\^E)dx (3.6) 

+ C j E\cuT\H\^dx. 

We estimate each term on the right-hand side of (j3.6p as follows: 
First, Holder's inequality gives 

j {\u\\Vu\+ pe\u\)\VE\dx <ri\\VE\\\2+C{ri) j [\u\'^\Vu\'^ + pE'^\u\'^) dx. (3.7) 

Next, if $7 = M^, Sobolev's inequality gives that there exists a universal constant C 
such that 

\\E\\l.<C\\VE\\l2. (3.8) 

If is a bounded smooth domain in M^, the Poincare-type inequality ( [S'l Lemma 3.2]) 
shows there exists a generic positive constant C which also depends on Vt such that 

II^IIlb < C\\p^/'^E\\l2 + C\\VE\\l2, 

which combined with ()3.8p implies 

||^||l6 < C\\p^''^E\\l2 + C\\VE\\l2. (3.9) 

It thus follows from Holder's inequality, ()3.9p . ()3.ip . and ()2.ip that 



j {\u\\H\'^\VE\ + \Vu\\H\^E)dx 



< Ml4h\\14^e\\l2 + \\Vu\\l2\\h\\14e\\l6 (^-^^^ 

< rj\\VE\\l,+C{rj)\\p'/'E\\l,+CmVu\\l,. 



Finally, integration by parts together with (13. 9p yields 

ElcmlHl^dx < C / \VE\\VH\\H\dx + C / |£;||V^i?||F|dx 



< C||V£;||i2||Vi?||i6||F||i3 + C\\E\\lc>\\V'^H\\l2\\H\\ls ^^-^^^ 

< r^\\VE\\l,+Cmp^/'E\\l,+CmVH\\l,. 
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Putting (13. 7|) . (I3.10p . and (13. lip into ()3.6p . we obtain ()3.4p after choosing ry suitably 
small. The proof of Lemma 13.21 is completed. 

Then, we derive the following crucial estimate on the L°°(0, T; L^)-norm of Vu. 
Lemma 3.3. Under the condition (|1.26p . it holds that for <T < T* , 

sup / {{p- pf + pO"^ + \Vu\'^ + \VH\'^)dx 

(3.12) 



Proof. First, multiplying (jl.ip 9 by ut and integrating the resulting equation over 
show that 

~ J {filVul"^ + {p + X){dwuf)dx + J p\u\'^dx 

= j pii- {u- V)udx + j Pdivutdx -\ j P - 2div(iJ H)) ■ utdx 



<^ j p\u\^dx + C j p\u\^\Vu\^dx + ^ j PdiYudx 
- j Ptdivudx + ^ j {\H\'^divut - 2H ■ Vut ■ H)dx. 



(3.13) 



Then, we will estimate the last two terms on the right-hand side of ()3.13p . 

On the one hand, to overcome the difficulty caused by the boundary, motivated 
by |12ll34j . we decompose the velocity into two parts. It follows from Lemma 12.31 that 
for any t € [0, T], there exists a unique v{t, •) G Dq n Z)^'^ n D^'^ satisfying 

pAv + {p + \)VdiYv = VP, (3.14) 

which together with (|2.5p yields that 

W^vWl, <C\\P\\lv <C\\pE\\l,, for [2,6], tG[0,r], (3.15) 

and that 



j Ptdivvdx = — J {pVvt ■ Vv + {p + X)divvtdivv)dx 
= -2jt j {l^\^v\'' + {p + \){diYvf)dx. 



(3.16) 



Denoting by 

w = u — V, (3-17) 
we have w G Dj^fl-D^'^nD^'^, for a.e. t £ [0,T]. Moreover, for a.e. t G [0,T], w satisfies 

pAw + {p + A)Vdiv?i; = pu + H x (cmlH), (3.18) 

which together with the standard L^-estimate for elliptic system gives 

llV-wll^e + ||V^i(;||i2 < C\\pu\\l2 + C|||F||ViJ|||£,2 

< C\\pu\\l2 + C\\H\\l4VH\\]^2\\VH\\]^^^ (3.19) 

< C\\pu\\l2 +C\\VH\\]I.^\\VH\\]^^, 
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due to (O). It follows from ([331) and SJM that 



— [ Ptdivwdx = —— [ {pE)tdivwdx -\ [ {p\u\'^)tdivwdx 

J Cv J 2Cv J 

<C j {pE\u\ + \VE\ + \u\\Vu\ + \u\\H\^) \V^w\dx 

f o Riy f o (3-20) 

+ C I \H\^\Vu\\Vw\dx - — I \cmlH\^dwwdx 

R I ^ 
— - — / (div(pii)|iipdiv?i; — 2pu ■ utdivw) dx = /«. 

Cauchy's and Sobolev's inequalities together with (13. ip yield that 

h+l2<v{\WM\l2 + \Ww\\l,) 

+ C{r]) j {p^E'^\u\^ + \VE\^ + \u\^\Vu\^ + \Vu\^)dx. ^^'^^^ 

Similar to (|3.1ip . integration by parts leads to 

h<C j {\VH\\H\\V'^w\ + \V^H\\H\\Vw\)dx 

<C [\\VH\\l4'^Ml^ + \W^H\\lA\'^M\l<^\\H\\l^ ^^-^^^ 
<ri + l|Vu;||ie) + CmVH\\\,, 

where in the last inequality we have used ()3.ip . 
Integration by parts also gives 

/4<Cy" p\uf\V^w\dx + C j {p\u\^\Vu\+ p\u\\u\)\Vw\dx 

<C{7]) j {pE^lul"^ + p\u\'^\Vu\'^ + p\u\'^\Vv\'^) dx (3.23) 
+ r/||V^?i;||^2 + J p\u\'^dx. 
On the other hand, direct calculations show 
j {\H\'^divut - 2H ■ Vut ■ H)dx 
= ^ j (li^l^divu -2H -Vu- H)dx 

-2 j \h ■ Htdivu - HfVu - H - H -Vu - Ht)dx 
- ^/ (l^l^divu-2F- Vn-F)(ia; + C||Fj||22 + C|||i/||Vn|||22. 



(3.24) 



Substituting (|3TB and ([3:20]) - (l3:2^ into (|3l^ . we obtain after using (IXT9]) and choos- 
ing r} suitably small that 

^ j ^dx + j p\u\'^dx 

<C j {pE^ + |Vn|2 + |Vz;|2) dx + C\\Vu\\l, (^"^^^ 
+ C2 {\\VE\\l, + WHtWl. + \\VH\\l. + |||//||Vn|||i.) , 
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where 



$ =/i|Vtip + (/u + X){divuf - 2Pdivu + /u|Vf p + + A)(divt))^ 



satisfies 



^ > ^|Vn|2-C3/9^2_c|F|l (3.26) 



Next, it follows from (jl.ll) that for r, s as in (11.7 

^|llv^lli2 + iii^*iii. + ^^llAi/|ii. 



(3.27) 



<C(|||F||Vn|||i. + ||H|Vi7|||i.) 

< r/llV^xllie + C(7?)||Vn||i. + C(e)(l + ||n||i,.)l|V//||i2 + ellV^i^H^,. 
Noticing that the standard L^-estimate of elliptic system gives 

llV^i^ll^, <C4||Ai/||^2, 
after choosing e suitably small, we deduce from (j3.27p that for any rj € (0, 1), 

4ujjVH\\l,+A\\H,\\l,+2u'C^' {\\VH\\l, + |||//||Vn|||i.) ^g) 
< Cr^WVuWl, + CmVuWl, + C{1 + \\u\\ir)\\VH\\l,. 

Then, adding 1^ multiplied by C5 = Csc^^ + (C2 + 2)k-^ and ([3:28]) by Ce = 



(1 + C4i^"^)(C2 + C1C5 + 2) to flOKI) . we obtain that 
^ y ($ + C5pS2 + 4C6Z^|Vi7|2)cix 

+ llVi^lli^ + p\u\'dx + WHtWl, + \\VH\\l, 

+ Cr^\\Vu\\l^ + Ci-q) j {pE^ + |Vn|2) dx. 
Holder's inequality together with (|3.15p yields that 

< C\\u\\lr {\\p'/'E\\l,^,,..,, + \\Vu\\ ^2r/(r-2) ) 

< C\\u\\l {\\p'/'Eft'''^\\E\\%^ + ||Vn||r^)/lVng^) 
<C(r?)(l + 11^x1110 y (pi^2^|V7x|2)dx + r/||Vi?||2,+,7||V^x||2„ 



(3.29) 



(3.30) 
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where in the last inequahty we have used (13.90 . It follows from ()3.17p . ()3.15p . p.l9p . 
and IK9\j that 



\\Vu\\l, < C\\pE\\Le+C\\pu\\r^+C\\VH\\Y^\\VH\\]^^, 

< C\\p^/^E\\l2 + C\\VE\\l2 + C\\pu\\l'^ + C\\VH\\li^\\VH\\^^l 



(3.31) 

L2 II \' Urn 



Putting this and ()3.30p into ()3.29p . and choosing rj suitably small, we obtain after using 
Gronwah's inequality, dSl]), and (fLMI) that 



sup / {pE"^ + |Vn|2 + |VF|2) dx 

0<t<T J 

+ j j (|V^p + p\u\^ + |n|2|Vnp + \Ht\^ + jV^Fp) dxdt < C. 



Finally, (jl.ip i implies that 

{p — p)t + div((p — p)u) + pdivn = 0. (3.33) 

Multiplying (|3.33p hy p — p and integrating the resulting equation over Q, we obtain 
after using ()1.26p that 

i\\p-p\\hnt)<c\\p-p\\i2+c\\vu\\i2, 

which together with (I3.32p and the following simple fact that 

\\^0\\l2 < C\\VE\\l2 + C|||n||Vn|||i2, 

directly gives (I3.12p . The proof of Lemma [3^ is completed. 

Finally, the following Lemma 13.41 will deal with the higher order estimates of the 
solutions which are needed to guarantee the extension of local strong solution to be a 
global one under the conditions (|l.lip - (11.13p and (ll.26p . 

Lemma 3.4. Under the condition (11.26j) . it holds that for < T < T* , 

sup {\\p - pWmnw^.i + + W'^SWm + WHWh^) < C. (3.34) 

0<t<T 

Proof. First, it fohows from (l33T]l . (fXT2]) . and 1^ that 

llV-ull^e < C + C\\pu\\l2 + C\\Ve\\L2 + C\\\u\\Vu\\\l2 + C\\VH\\]i^ 

<C + C\\pu\\l^ + C\\Ve\\L2 + C\\u\\Le ||Vu||J^/^||Vn||^/^ + C\\VH\\]i^ 



<C + C\\pu\\l2 + C\\Ve\\L2 + ^||Vn||i6 + C\\VH\\]il 



which implies 



Vn||i6 <C + C\\pu\\l^ + C\\Ve\\L2 + C\\VH\\^fil (3.35) 



I// 

Then, it follows from the standard L^-estimate of (jLip ^. (|3.12p . and ()3.ip that 

\\V^H\\l2 < C\\Ht\\L2 + C\\\u\\VH\\\l2 + C\\\Vu\\H\\\l2 

< C\\Ht\\L2 +C\\u\\L4'^H\\\f,^\\V^H\\]^,^ + C\\H\\l4\\Vu\\l4 

< C\\Ht\\L2 + ^\\V^H\\l2 + C||V^x|U4 + C, 
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which together with (I3.12p and ()3.1|) imphes 

\\H\\h2 < C\\Ht\\L2 + C||Vn||i4 + C. (3.36) 
Holder's inequahty, along with (I3.35P and ()3.12p . gives 

||Vu||i4 < ||Vn||J^//||Vu||5^// < C + C\\pu\\%^ + C\\V9f/2 +C\\VHfj^t (3.37) 
which combined with (I3.36P and (I3.35P shows 

||Vn||i6 + \\H\\h2 < C\\pu\\l2 + C\\Ht\\L2 + C||V0||i2 + C. (3.38) 
Then, similar to ()3.9p . we have 

||^||l6 < c\\p^/^9\\l2 + c\\ve\\L2 <c + c\\ve\\L2, (3.39) 

which together with the standard L^-estimate of (jLip ^j and p.l2p gives 

\m\'m 

<C + C\\Ve\\l2 + C j pe^dx + C j p^e^\Vu\'^dx + C\\Vu\\l^ + C\\VH\\\i 
<C + C\\Ve\\l2 + C j pfe + C7||Vn||^2||0||2^ +C||Vu||^4 + C||VF||^i 
<C + C j pe^dx + CllV^ll^a + ^llV^ll^i + C||Vn||^4 + C\\VH\\%,. 
Combining this with (|3.36p shows 

llV^ll^i <C J pe^dx + C\\Ve\\l2 + C||Vn||^4 + C\\Ht\\l2 + C. (3.40) 

Next, we claim that we have the following estimates on both u and 9, (j3.4ip and 
()3.42p . whose proofs are similar to those in [15''38] and can be found in Appendix A 



sup [ (\Ve\^ + p\u\^ + \Ht\^)dx+ [ [ (pe^ + \Vu\^ + \VHt\^)dxdt<C, {3A1) 

<t<T J Jo J ^ ' 

sup ||p^/^^||i2+ / \\Vef^2dt<C. (3.42) 

<t<T Jo 



0<t<T 

Then, the combination of (fOSjl - ffS:!?]) with (fS:^ leads to 



sup {\\Vu\\Le + \\H\\H2 + \\e\\L6 + \\Ve\\Hi) <C. (3.43) 

0<t<T 

For 2 < p < q, direct calculations show that 

^IIVpIIlp <C(1 + \\Vu\\lo.)\\Vp\\lp + C\\VMlp. (3.44) 
For V £ D^n D"^'^ D D^''' satisfying (fXTip . it fohows from Lemma [231 and (pHHj) that 

\\Vv\\l^ <C{1 + log (e + ||V(p^)||^,-) IIp^IIl- + \\pe\\L2) 

<C71og(e + ||Vp||^,-). ■ ^ 
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Then, foi w = u — v £ Dq H D^'^ n D^'^ satisfying (I3.18p . applying the standard 
LP-estimate to ()3.18p . along with ()3.43p . gives 

\\V^w\\l6 < C||pn||i6 +C7|||i/||V//|||i6 
< C\\Vu\\l^+C, 

which together with ([XT9D . ([3l3t\ . and ([OTD shows 

||Vi(;||loo < C + C||Vtt||i2. 

The combination of this with (I3.45P gives 

IIVuIIloo < Clog (e + ||Vp||i,-) + C\\Vu\\l2. (3.46) 

Applying the standard L^-estimate to ()1.1|) 9 leads to 

||V\||lp < C(||pu||lp + |||i/||Vi/|||LP + ||VP||lp) 

<C{\\pu\\Lp + \\yp\\Lp) + C (3.47) 
<C{l + \\Vuh^ + \\Vp\\Lp), 

due to and (ITIHD . Substituting (ITTTD and into (IXH) yields that 

f\t)<Cg{t)f{t)lnfit), (3.48) 

where 

/(t)^e+||Vp||^,-, 5(t) = l + l|Vii||L2. 
It thus follows from ()3.48p . ()3.4ip . and Gronwall's inequality that 

sup \\Vp\\l,<C, (3.49) 

0<i<T 

which, along with p.46p and (13.4ip . directly gives 



T 



\Vu\\lo.dt < C. (3.50) 







Taking p = 2 in <^JM^ we get by using ([330]) . ([3311), (fOT]) . and Gronwall's inequality 
that 

sup \\Vp\\l2<C, (3.51) 

0<t<T 

which together with (l3:i7D . (ITiH]) . and (IHIHD yields that 

sup \\V^u\\l2<C sup (||/9tt||i2 + ||V/9||i2 + ||V0||i2 + ||// • V-H'||i2) < C. 

0<t<T 0<t<T 

This combined with ([OO]) . ([33T]) . (jOS]) . and ([3T2|) finishes the proof of Lemma [331 
Now we are in a position to prove Theorem ll.il 

Proof of Theorem ] Suppose that (I1.14p were false, that is, (I1.26P holds. Note that 
the generic constant C in Lemma 13.41 remains uniformly bounded for all T < T*, so 
the functions {p,u,9, H){x,T*) = lim {p,u,6, H){x,t) satisfy the conditions imposed 
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on the initial data (jl.lip at the time t = T* . Furthermore, standard arguments yield 
that pii, p9 G C([0,T];L^), which implies 



{pu,pe){x,T*) = lim^{pu,p0) e L . 



Hence, 



-fiAu -{p + A)Vdiv'u + RVipO) - (curl H) x H\t=T* = ^/p{x,T*)gi{x) 
kA9 + ||Vu + (Vn)*n^ + X{dwuf + v\cuvl H\^\t=T* = ^{x,T*)g2{x), 



with 
and 

g2{x) -- 



A 



' p-^/'^{x,T*){pu){x,T*), for X e {x\p{x,T*) > 0}, 
0, for X £ {x\p{x,T*) = 0}, 



p-'^/'^{x,T*){c^pe + Rpedwu){x,T*), for x £ {x\p{x,T*) > 0}, 
0, for X G {x\p{x,T*) = 0}, 



satisfying gi,g2 G due to (^Mi), (1332]), and (|33i]l . Thus, {p,u,e, H){x,T*) also 
satisfies ()1.12p and (jl.lSp . Therefore, one can take (p,u,9, H){x,T*) as the initial data 
and apply Lemma l2.1l to extend the local strong solution beyond T* . This contradicts 
the assumption on T* . We thus finish the proof of Theorem 11.11 



Appendix A Proofs of and ( KW i. 



The proofs of ()3.4ip and ()3.42p are a direct combination of those of Lemma 4.1 and 
(4.28) in [15] with that of (3.24) in [38]. We sketch them here for completeness. 

First, it follows from (f3T2]l and (f3?38]) that 



sup / p9^dx + / {\\V9\\l2 + llVnll^e) dt < C. (A.l) 

^<t<T J Jo 



2di 



_ / o9^dx + / 

0<t<T. 

Applying Uj[dt + div(ii-)] to (jl.iP n and integrating the resulting equality over give 
J p\u\'^dx = — J Uj[djPt + div{udjP)]dx + 1^1 J Uj[dtAuj + div{uAuj)]dx 
+ (^ + A) J Uj[djd\vut + d\v{udjdivu)]dx 

iijldtdjlHl"^ + div{udj\H\'^)]dx 2) 
+ J Uj[dtdi{WHj) + d\w{udi{H'Hj))]dx 

5 



i=l 



19 



We get after integration by parts and using the equation (jl.ip i that 
Ni = - J Uj[djPt + dW{d jPu)]dx 

= i? y djiij (^p6 — pu ■ V6 — 6u ■ V p — 6pd\vuj dx + j dkUjdjPukdx 

Opdivuj dx + J Pdiviidivudx — J PdkUjdjUkdx ^ 



R / d. 



< f llVn 
<^l|Vn 



I2 + c\\pe\\l, + c J p^e^\vu\^dx 
i. + c||p^|ii. + c||/>0||^/'ll^ll^//livu||i4 



\h + c\\pnh + c\\ve\\l, + c\\vu\\l, + c, 

where in the last inequahty we have used (|3.39p . Integration by parts leads to 

N2 = p j Uj[dtl^Uj + di\'{ul\Uj)\dx 

= —p J {diUj{diUj)t + AujU • Viij) dx 

= —p J [\Vu\'^ — diUjUkdkdiUj — diiijdiUkdkUj + AujU ■ Viij^ dx (A. 4) 

= —p J (iV-itp + diiijdiUjdivu — diiijdiUkdkUj — diUjdiUkdkUj) dx 
<-^ / 1^^!^^^ + '^ / \Vu\^dx. 
Similarly, we have 

A^3 < -\{p + A)||divu||i2 +C j \Vu\^dx. (A.5) 

Integration by parts together with (j3.12p and (|3.ip shows 

m < C\\Vu\\L2{\\\H\\Ht\\\L2 + \\\u\\H\\VH\\\l2) 

< C\\Vu\\L<\\H\\L4Ht\\]!2\\Ht\\]!,^ + \\u\\l4H\\l4'^H\\l'^) (A-6) 

< e||Vn||i. + vW^mh + C{e,v)\\Ht\\l, + C{e)\\V^H\\l,. 

Similarly, we also have 

liVsl < e\\Vu\\l, + r?||Vi/,||i. + C{e,v)\\Ht\\h + C{e)\\V'H\\l. (A.7) 
Substituting (IXSjl - dATT]) into ^A2\i . we obtain after choosing e suitably small that 

1 1 p\n\'dx + p\\Vu\\l,<C I p9^dx + Cv\\VHt\\l,+CmHt\\l, 

+ c\\v^H\\l2 + c\\ve\\l2 + C\\Vu\\l, + C. 
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Next, multiplying (ILlh ci by 9 and integrating the resulting equality over yield that 

^{\m\h)t + cvl P\0\'dx 
= -K jve- V{u ■ Ve)dx + \ j {d\wufedx 

+ \^{u)\'^edx - R j pOdwuMx + v j \cm\H\'^edx 



(A.9) 



5 

A 



i=l 



We estimate each = 1, • • • , 5) as follows: 
First, it follows from dOOjl and (f3l^ that 



|/i| <c j \Vu\\ve\^dx 



<c\\vuU4ve\\'/^\\ve\tl^ 
<5\\vH\\l, + cm'^e\\l, 



(A.IO) 



<C5 j pO^dx + C{5)\\Ve\\l2 + C||Vn||^4 + C\\Ht\\]^^ + C. 

Next, integration by parts yields that, for any rj E (0, 1], 
/2 =A j {d\Nufetdx + \ j {d\wufu ■ V9dx 

=X j {{dwufe)^ dx-2X j edivudiviii - u ■ Vu)dx + X j {divufu ■ VOdx 
=A {dwu)'^9dx^ — 2X J Odivudivudx 

+ 2X J OdiwudiUjdjUidx + X j u-V (6'(divu)^) dx 

<x {^j {d\Yufed3^ + c\\d\\LA\yu\\^ff\\^utl! (||Vu||i2 + \\vu\\l,) 

<X {^j {dbjufOdj!^ + r/llVnlli^ + C(??)||Vn||^4 + C||V0||^2 + C7, 



(A.ll) 



where in the last inequality we have used ()3.39p . 

Then, similar to (jA.lip . we have that, for any r/ G (0, 1], 



l3<2p(^j \^iu)\^9dx^ +r,\\Vu\\l2 + C{r,)\\Vu\\l^ + C\\V9\\l2 + C. (A.12) 



Next, it follows from (l3J2]l and dOOjl that 

< 5 J pe^dx + Ci5)\\Ve\\l2 + C7||Vn||^4 + C{d), 



(A.13) 
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and that 



h = T^ j \cnT\H\^etdx + v j \cm\H\'^u-Vedx 

= ^-^j \cnT\H\'^edx -2v j OqmiIH ■ cnilHtdx + u j \c\it\H\'^u -S/Odx 

+ c\\VH\\l4u\\LA\^eu. 

<ujj |curlF|20dx + vWVHtWl, + C{rj) (l + ||V0||i.) (l + \\V^H\\l,) . 

Substituting ([O0]) - (|Al4l) into (jX9]l . we obtain after choosing 6 suitably small that, 
for any r] £ (0, 1], 



"^dx + Ou I p\9\'^dx 



d 

di _ 

< C7(r?) (1 + llV^lli.) (1 + \\V'H\\l, + II V^lliO + CvWVuWl. ^^'^^^ 
+ Cr^WVHtWl^ + C||Vn||i4 + C\\Ht\\l, + C, 



where 

4 k|V0|2 - 29 [A(divn)2 + 2^|D(n)|2 + z^|curlFp] . (A.16) 

Next, differentiating (jl.ip /i with respect to t and multiplying the resulting equations 
by Hf, we obtain after integration by parts and using (13.ip and (13.12P that 



j \Ht\^dx + v j \VHt\'^dx 



1 d 

<C(|||nt||/7|||i2 + |||n||/7t|||i2)||V//j||i2 

< C (llliillii'lllLa + II |^/||V7x| 1^111^2 + \\u\\L4Ht\\]!2\\Ht\\]!,^^ W^Hth^ 

< C (||n||,.6 11^11^3 + 11^x11^6 II V7x||i4||//||ii2 + ||i7t||^/'||V/7t||^/') \\^Hth2 



< 



V 



2 

which implies 



||V//t||i2 + CllVtilli^ + C||Vtx||i4 + C||Fi||i„ 



^ y \Ht\^dx + iy J \VHt\^dx < C||Vii||i2 + C\\Vu\\l, + C\\Ht\\l2. (A.17) 

Finally, adding (|A.8j) multiphed by r/^/^ and (lA.17j) by t?^/^ to (lA.lSj) . we obtain 
after choosing t] suitably small and using (|3.37p that 

2^ J\^+v'/W+v'^'p\n\')dx 

+ J (^o,p\9\^ + u7]^/^\VHt\^ + fir]^/^\Vu\^'^ dx (^.18) 
<C(r?) (l + ||V0||i.) (1 + ||V2/?||2, + ||V0||2^) 
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Noticing that (lAl6]l . (l3J2]l . (IM]) . and (I09D lead to 

> 2A.||v0||i. - c\\eu4vu\\]^^\\vu\\'/^ - cpu.wvHW^ffwvHW^J^l 

+ 2 l{rj'/^\Ht\' + ri'/^p\u\')dx 

> K||V0||i. + y" (r?i/2|//^|2 + riy^p\nf)dx - Civ), 

we directly obtain p.4ip after using Gronwall's inequality, ()1.12p . (jA.lSp . (|3.12p . and 
It remains to prove ([232]) • First, it follows from ([335]) -([331]) that 

sup (11^11^6 + ||Vn||i2nL6 + \\H\\h2) + T \\V^e\\l,dt < C. (A.19) 

0<t<T JO 

Similar to ()3.9p . we have 

\\e\\L'^<c\\p''H\\L^ + c\\ve\\L2. (A.20) 

Next, applying the operator dt + div(ti-) to p.ip -^ leads to 
c^p (dt6 + U-V6 



= kA0 + K (divnA6' - dt [diU ■ VO) - ^^u ■ VdiO) 

+ (A(divn)^ + 2^|D(n)p) divn + RpddkUidiUk 21) 
— RpOdivu — RpOdvvu + 2 A (divii — dkUidiu^) divu 
+ p{diUj + djUi) {diiij + djiii - diU^dkUj - djUkdkUi) 
+ u [at|curlF|2 + div(|curl/7|2n)] . 

Multiplying ()A.2ip by 9, we obtain after integration by parts and using ()A.19p . ()3.4ip . 
and ([A?20]) that 

'-^ 1^1 p\e\^dx^ + n\m\h 

<C j |Vn| (iV^^II^I + |V^||V^|) (ij; + C7 j \V u\'^\e\ {\V u\ + 6) dx 
+ C j p\9\'^\Vu\dx + C j pe\Vu\\e\dx + C j \Vu\\Vu\\e\dx 
+ C j {\VH\\yHt\\6\ + \VH\'^\u\\Ve\^dx 

< c\\VuU^\\ve\\H^ (\\e\\L. + WvHl^) + c\\Vu\\l4e\\L. {\\VuU^ + pu,) 
+ c\\Vu\\M\LAmL^+c\\p''H\\'^^^^^ 

+ c\\Vu\\LAWu\\L4e\\Le + c\\VH\\L4VH4ML^ 
+ c\\VH\\l4u\\L4'^eU. 

< ^\m\l, + C\\V'e\\l, + C\\p'/'9\\l, + C||Vn||i. + CWVHtWl. + c, 

which combined with Gronwall's inequality, (jl.lSp . ()A.19p . and (|3.4ip directly gives 
(1332]). 
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